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Summary 

The dynamic motion of a cylinder with and without end caps is analyzed based on rotation about its 
center of mass and restrained by a plane normal to the axis passing through its center of mass at an angle 
a. For small values of a, the governing equations are simplified, and for symmetric bodies, stability 
requires rotation greater than [2 ^l(JWL*)\/J x , where J is the transverse mass moment of inertia, W is the 
weight of the cylinder, L* is the cylinder length from the base to the center of mass, and J x is the mass 
moment of inertia about the longitudinal axis OX of the barrel. Comparisons to data are made and some 
applications are discussed. 


Introduction 

The dynamic motion of a cylinder on a floor or hard plane surface is entertaining and instructive. 
Many have witnessed the rolling and collapsing motions of a runaway tire, yet the similar behavior of an 
oil drum is rarely seen. The dynamics of the latter are difficult to describe as the motion can be induced 
into the drum by a regular periodic change in torque on the periphery of the drum, such as once every 
time the drum rotates about its transverse axis. The major dilemma resides in its control, which seems to 
be an accomplishment of those skilled in rhythm and agility. Even though difficult to describe and put 
into practice, such motions are of interest to rotor dynamicists in stabilizing and maintaining stable 
spinning spacecraft and in spontaneous unloading of transmissions. 

Although the outcome compares favorably with known visual data, we recognize several limitations 
in that maintenance torques, friction, and other dissipation mechanisms are not completely addressed. The 
current analysis is not dedicated to the transient period (shown in the video clip which can be viewed in 
the PDF file of this document at http://gltrs.grc.nasa.gov/reports/2006/TM-2006-213583/TM-2006- 
213583.pdf). We are mostly concerned with the stability of the barrel after a stationary situation has been 
reached. Further, while ignored herein, the issues of maintaining control of a large rotating mass are 
clearly the key, yet can be effectively accomplished with minor periodic adjustments. Flowever, the 
present report contains only a preliminary analysis of these types of motions. The objective of this 
analysis is to obtain analytical stability conditions; thus, the model contains inherent approximations and 
simplifications. Future reports will address different aspects that are not treated herein. 

The appendix contains a summary of the formulas that were used to calculate the inertial 
characteristics of the drum. 

Figure 10 illustrates the dynamics and control of a rotating barrel that is set in motion, lifted, and 
stabilized in spinning about the longitudinal axis via a series of fast and relatively weak impulses applied 
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by a skilled operator. A video of this motion can be viewed in the PDF file of this document at 
http://gltrs.grc.nasa.gov 


Symbols 

A parameter defined in eq. (28) 

A * parameter dependent upon sliding conditions 

a constant coefficient in eq. (42) 

C\ C 2 generic notation for integration constants to be determined from initial conditions 
c constant in eq. (42) 

D diameter, m 

Ff friction force 

H height of barrel, m 

J transverse mass moment of inertia, kg m 2 

Jx mass moment of inertia about longitudinal axis OX of barrel, kg m 2 

Kq angular momentum, kg m 2 /sec 

L length of barrel, m 

L* length of cylinder from base to center of mass, m 
M total external moment, Nm 

m mass of disk 

N normal force 

R radius of barrel, m 

r position vector of the acting point of a force, m 
t time, sec 

W weight of cylinder, N 

a nutation angle, deg 

q constant in eq. (44) 
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0 polar angle in a set of cylindrical coordinates (r, 0, z ) 
p friction coefficient 

v precession angle, deg 

^ location of center of mass measured from base upward in fractions (0 to 1) of half of the 
length of the barrel 

p density, kg/m 3 

a constant parameter in eq. (42) 

a s surface element in eq. (56) 

x thickness of metal sheet, m 

cp Eulerian rotation 

X geometric parameter ( = L * sin a - R cos a ), m 

'F dummy parameter, where 'Ey = 'F along /-coordinate axis, 'F t = tbF/ck, *F rt = d 2v F/dr 

co angular speed, rad/sec 

Subscripts: 

CG center of gravity 

N normal force 

p precession 

S in R s represents radius of circle described by the base of the spinning drum 
t time derivative 

The main geometric elements are indicated in figure 1. 

The position of the mass center is L * = qL/2, where \ is the location of the center of mass measured 
from the base upward (in percentages of half the length L of the barrel). 
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Figure 1 . — Position angles and rotations of frames of axes. 


Frictionless Analysis 

Preliminary Aspects 

The dynamics of the cylinder is the same whether the cylinder has closed ends or not. With fixed ends 
or heads, the cylinder is termed a “drum,” without a closed upper end, an “open drum,” and without either 
end closed, a “plain cylinder.” In all cases, the restraint is that of a flat surface upon which the cylinder is 
supported and rotates. To support the motion, the cylinder must spin on multiple axes, and this is the 
motion analyzed herein. 

The dynamic modeling of a body can be done in various ways. The method and the position variables 
can be selected as desired. The analysis developed herein is governed by the following assumptions: 

(1) The drum is modeled as a symmetric rigid body whose center of mass is fixed during the motion 
and is restrained from below by a plane surface. 

(2) The center of mass is located on the centerline of the barrel. 

(3) The motion is described by the theorem of the moment of momentum or the angular moment. The 
equations will be written in the body-fixed frame OXYZ. 

(4) The frame OX F Y F Z F is fixed in space. 

(5) The frame OXYZ is fixed with respect to the drum and executes all the motions of that drum. 

The order of rotations necessary to obtain the OXYZ frame from OX F Y F Z F is (fig. 1) 

v about OX F = OX\ 

a about OZ\ = OZ 2 

cp about OX, = OX 
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which is sometimes schematically written as 


OX F Y F Z F 


OX, =ox. 


o()X x Y x Z v 


a 


OZ| =OZ 2 


xOX 2 Y 2 Z 2 


9 


OX 1 =ox 


-xOXYZ 


The friction is not taken into account in the initial analysis; however, it will be added at a later stage. 
When the friction is neglected, the only moments that act upon the body are produced by the normal force 
N exerted by the plane. 

The components of the normal force are (fig.2): 

N x = W cos a 

Ny — ~W sin a cos cp (1) 

N z = W sin a sin cp 

The components of the angular speed along the body-fixed axes OXYZ are (fig.3) 


<x>x = + v t cos a 

G) Y = a t sin cp - v t sin a cos cp (2) 

co ^ =u t cos cp + v^sina sin cp 

and the components of the position vector of the acting point of the normal reaction force r are (figs. 1 
and 2) 



Figure 2. — Force components and position vectors. Friction neglected 
for case 1 . 
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Figure 3. — Components of angular speed. 

{ r Y) N = Rcos( P (3) 

{ r z) N 


Governing Equations 

The moment-of-momentum equations around the fixed-mass center are written in terms of the body 
axes OXYZ : 


d K, 


o 


/■ =J — - + coxA: 0 l = {m} 

d t dt 1 

J XYZ l J XYZ 


(4) 


When the friction is neglected, the only moments are produced by the normal reaction force 


M -r N xN = 


i j k 

-L* .R cosrp -.Rsinrp 

W cosrp -W sin a cos rp W sin a sin rp 


(5) 


so 


M x = 0 

M Y = ifsinrp(L* sina-i?cosa) = JTxsinrp (6) 

M z = W cosrp(L* sin a -i?cosa) = W% cosrp 
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Next, 


and 


^0 - J(0 Y 




0 

< Jx^x^Z ~ 

TcOj^COy — J X^X^Y 


Collecting the terms, the OX equation is 


Jx^Xt ~ M X 


So 


and further 


cojf = constant 


cp^ = U) X - v t cos a 


The OY and OZ equations are 


Jo) Yt +{J X -M Y 

Jo) Z t +{J~ Jx) (0 X (0 Y = M Z 


Next, 


co Yt = ct tt sincp + oqcp ? coscp- v tt sin a cos cp- v t a t cosacoscp + v ? cp ? sin a sin cp 


co z t = a tt cos <p-a f <p f sin cp + v tt sin a sin cp + v t a t cos a sin cp + v t q> t sin a cos cp 

Next, equation (13) is multiplied by (-cos cp), equation (14) is multiplied by (sin cp), and the 
expressions are added. Using equations (12) results in 

My J y ~ J 

-oqcp ? + v tt sin a + v t a t cos a = — coscp + — a>jU°z coscp 


M 7 . J — J y 

+ —j- smip j co x coy sin cp 


Using equations (6) and then equations (2) yields 


J X ~ J 

-a ? cp ? + v tt sma + v ( a ( cos a = — co^a t 


(?) 

( 8 ) 

(9) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 

(16) 
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For v, = constant, that is, v n = 0, the above expression becomes 


-a ? cp ? + v t a t cos a : 


l J 


® x a t 


(17) 


Now, the body in stable dynamic equilibrium executes stable oscillations about the equilibrium 
position. So a, is non-zero and it simplifies. Also, cp, from (11) yields 


v,= Jx ™ x 
2Jcosa 


(18) 


and further 


<Pt 


f Jx A 


1 -- 


V 


2 J 


(19) 


Next, equation (13) is multiplied by (sin cp), equation (14) is multiplied by (cos cp), and the two 
expressions are added. Further, using equations (12), then equations (6), and finally equations (2) yields 


W / ± ■ \ J Y 

a H +v ? cp ? sma = — (Z sma-iccosa)- 


J y — J 


J 


-co x v t sina 


( 20 ) 


From cp, of (1 1), the above expression becomes 


W , 


J i 


2 - yv ( t*. - \ J y 

a tt — v t sin a cos a = —[L sma-ACOsaj-— y-co^v^ sma 


( 21 ) 


and finally, with v, of (18), the governing equation of the a angle is 

W 


j 2 ^ 

a tt + — ^-co^ tana = sina- A cos a) 


( 22 ) 


For small a, tan a = sin a = a and cos a = 1; thus, the equation becomes a constant coefficients 
equation: 


( r2 


«« + 


Jx 2 W T * 


4 J- 


-c°^- y Z 


a- — 


WR 

J 


(23) 


The solution a = (Xhomogeneous + ^particular is stable when the homogeneous part does not contain 
exponentials of the form exp(±/f), with r a real number. So, we need the characteristic equation 


r 2 + 


( t 2 


J X 2 W * 

—UX-—L 


4 J 


J 


= 0 


(24) 


to have complex solutions such that 
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which yields the stability condition 


2yJWL* J 
co x > 

J X 


The analytical solution of equation (23) (for small a angles) can be written as 


WR _ . 

a = + a sm 

JA 



(26) 


(27) 


where 


A = 


^ Jx c °x 

4 J 2 


WL* 

J 


- 

( WR) 

, a = J 

a o + ~rr 

J v 

v JA J 


an, . _1 

H , and cp = tan 

A 


a 0 + 


WR 

~JA 


iVT 

a 0l 


(28) 


It can be readily seen that for zero initial conditions (a 0 = 0 and d 0 = a ()/ = 0 ) and for W —> 0, the motion 
is not an oscillation and the axis of the body remains oriented toward a fixed position. This is in fact a 
well-known property of a gyroscope, which is not subject to external influences (also known as Euler 
gyroscope). 


Simplified Analyses 

In an alternate, much simpler perspective, just by analyzing the static neutral stability with no slip 
interface for a+ = arctan(2/(77?) and Rs = R sec a+, if 

n d(p . , _ dv . .. dv .... 

R — = nm speed = Ac — and it — = constant (29) 

dt * dt dt 

and we have pure rolling with no slip, then if a increases, dip/ dr and Rs increase. Also, if a decreases, both 
dtp/dr and R s decrease. 

Also, in an even more simplified analysis, the balance point on the rim of an ideal cylindrical drum, 
albeit unstable, is at an angle 


a ^ = tan 


' 2R ' 

Ji, 


with an equilibrium rolling radius of 


R s -Rseca/, 


(30) 


( 31 ) 


For a barrel, one may approximate that H = 0.87 m and 2 R = 0.58 m. So, a* is about 0.6 rad and i?sis 
about 0.8 R. More accurately, a b = 0.588 rad (33.7° ) and R s = 0.2413 m 
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Alternate Approach 


No close solution can be found for the general equation of the a angle, and a complete study of the 
stability can only be done numerically. Many bodies with a gyroscopic motion are generally known to 
become unstable at high nutation angles, so the small-angle stability condition is usually a good indicator 
for the stability domain. If a more detailed stability analysis is desired, equation (22) has to be integrated 
numerically. Figures 4(a) and (b) compare the analytical solution for small a with the numerical solution 
for a closed drum with a L/D ratio =1.5 for cox three and four times higher than the theoretical stability 
limit, respectively. Figures 4(c) and (d) are similar to figures 4(a) and (b), but they are obtained using a 
slightly larger initial angle oto. Now, the current analysis is not aimed at obtaining a complete equation for 
the nutation angle. As it was already mentioned in the Introduction, herein we are only concerned with 
stability. The current modeling is not intended to provide a complete equation that would determine the 
accurate value of the angle a. The equation obtained herein can only be used as an indication of the 
stability of the motion. A stable solution of the equation can just indicate a stable behavior of the body, 
whereas an unstable solution is a warning that the body is unstable. 

For the completeness of this report, it is worth recalling that the analytical analysis presented above 
contains inherent approximations and that an accurate study of the dynamic requires the integration of the 
full version of the governing equations. The computation of the actual values of the variables requires the 
use of the angular momentum equations (12) along with the kinematic equations that can be readily 
obtained by solving equations (2) for cp, ,a„ and vy 

Figures 4(e) and (f) present results obtained by the integration of the general equations. Here the full 
system of the governing equations is integrated numerically. The variable cry is gradually increased until 
it reaches a steady-state value equal to three times the critical value determined by equation (26). The 
nutation angle oscillates steadily around a median position, as shown in figure 4(e). The orbit of the 
center point of the top lid is shown in figure 4(f). This orbit shows the series of fast small oscillations that 
are “wrapped” around a circle. To conclude, the analytical model is a good indicator of the stability; 
however, a complete dynamic analysis requires the full set of equations. 

Equation (22) can be integrated using the standard numerical procedures, such as Runge-Kutta or 
Adams. An alternate solution method can be also considered. 

First, let a, = C, , and then equation (22) becomes 


w 

da 


+ a\ tan a = a 2 sin a - a 3 cos a 


where a x = [( j\ /(4J 2 )] 0 ) 2 X , a 2 = ( W/J) (£, L/2 ), and a 3 = ( WtJ)R . 
Next, 


so 


and 


1 2 

—C,-a\ In (cos a ) = -a 2 cos a - a 3 sin a + Cj 


= [<q In (cos a) - a 2 cos a - a 2 sin a + Q 2 


(32) 


(33) 


(34) 
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Nutation angle, a, rad Nutation angle, a, deg Nutation angle, a, deg 



0.0 0.5 1.0 1.5 0.0 0.5 1.0 1.5 

Time, sec Time, sec 



Time, sec 


Figure 4. — Analytical and numerical solutions, (a) Speed, m x - 3*[2 'l(JWL*)]IJ x . (b) Speed, m x - 4*[2'I(JWL*)]IJ X . 
(c) Speed, co x = 3*[2^I(JWL*)]IJ X . (d) Speed, co x = 4*[2^ (JWL*)\IJ X . (e) Complete model integration; speed, 
(co^steady = 3*[2 ^l(JWL*)]IJ x . (f) Complete model integration; orbit of center of top lid (top view); speed, 
(^steady = 3*[2 V(JI/ML*)]/J X . 
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(35) 


t — C r > + 


|[«1 In (cos a) - a? cos o' -t /3 sin a + C] ] da 


Next, let £ = cos a to yield 

f = C 2 -J C 1 +a 1 lnC-a 2 C-a 3 (l-C 2 ) 


1/2 


n-1/2 


- 1/2 


1-C d£ 


(36) 


Again, the expressions have to be computed numerically. 


Applications (Computations for Barrels) 

Several applications are presented for the cases of an open barrel, a completely closed barrel, and a 
barrel with the bottom lid closed and the top lid removed. The typical barrel was 0.58 m in diameter, 
0.87 m in height, and had a metal thickness of 1.01 mm. The density of the steel was 7753 kg/m 3 , so the 
weight of the closed drum was 163.3 N. For this drum, the LID ratio was 1.5. The mass and inertia 
moments are presented in the appendix. We considered other cases for L/D ratios of 1 and 2. The results 
are summarized in table 1 . 


TABLE I.— THEORETICAL STABILITY LIMITS 


[Diameter, D, 0.58 m; minimum stable speed, coy, rad/sec.] 



For length-diameter ratios LID of — 

1.0 

1.5 

2.0 

a>x, rad/sec 

W, N a 

coy, rad/sec 

W, N a 

cojy, rad/sec 

W, N a 

Drum 

1.38 

122.5 

20.3 

163.3 

27.8 

204.1 

Open drum 

10.8 

102 

16.7 

142.9 

23.5 

183.7 

Plain cylinder (shell, no ends) 

10.6 

81.6 

15.9 

122.5 

22.3 

163.3 

Static angle a h , deg 

45 

45 

33.7 

33.7 

26.6 

26.6 

Equilibrium R s , m 

0.41 

0.41 

0.35 

0.35 

0.32 

0.32 


T Newton = 0.224808943 lb f 


Effect of Friction 

Two types of friction act upon the drum: rolling friction is perpendicular to the OY 2 axes and opposes 
the precession of the drum; sliding friction is parallel to the OY 2 axes and opposes the increase of the 
nutation angle a (fig.5). The rolling friction is much smaller than the sliding friction, so in the subsequent 
analysis is neglected (in industrial gyroscopic applications this is usually compensated by an electric 
motor). 

The components of the sliding friction force are 

F jx - Fj- sin a 

Fjy - iy cos a cos cp (37) 

Fj 2 = -Fj- cos a sin cp 

and the friction force is applied in the same point as the normal N. Following the same development as 
before, the equation for the angle a becomes 
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X2=X 4 



Figure 5. — Forced components. Sliding friction included for case 2. 



In a general case, the friction force can be written as 

iy = | jlN + A*a t 

where |i is the coefficient of friction and A * depends upon the sliding conditions. So 



(38) 


(39) 


(40) 


(41) 


where p is the friction coefficient. 

Now, the friction force in the left-hand-side term is usually much smaller than the other terms, so it 
will be neglected. As a consequence, the stability condition does not change and the equation can be 
written as 

9 

a tt + aa t + co“ aa = c (42) 
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where the constants a, c, a, and o) /; can be readily identified by inspection. For example, 


co 


2 

P 



J X (a X 
4 J 2 


and ct = 1 - 


WL* 

Jafp 


(43) 


The solution is 


sin(c o p ^at + C 2 ) + -^— (44) 

CO pO 

where C\ and Ci are integration constants to be determined from the initial conditions and p = 
a 2 /^4co^aj-l. 

It can now be readily seen that for a > 0, the oscillatory component decreases rapidly and the solution 
may appear to be a constant equal to the particular solution of the nonhomogeneous equation; that is, 

a = (45) 

CO "pO 

However, the result provided by equation (27) may not be the correct steady-state value. The correct 
steady-state value can only be obtained by integrating equations (2), (4), and (5). 


a = Ci 


exp 


^ -at ^ 


Additional Applications 

The foregoing treatment of the barrel follows general methodologies that are usually applied for 
gyroscopic motions. Other different applications can be found in various sources. Maybe the most simple 
and the most intuitive is the steady precession of a symmetric gyroscope or top (fig. 6). The motion is 
described as shown in the Governing Equations section; however, here the angular momentum theorem is 
written with respect to the contact point. Again, in a first approximation, the effects of the friction may be 
neglected, so the only moments that act upon the body are produced by the weight. 

The components of the load are 


W x = -W cos a 
W Y -W sin a cos cp 
W z = -W sin a sin cp 


and the coordinates of the mass center are 


(fy) w =L* sin a cos cp 
( r z ) w =-L* sin a sin cp 


(46) 


(47) 


So again, 
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Figure 6. — Steady precession of symmetric top. 


M x = 0 

M Y -WL* sin a (1 + cos a) sirup 
M z = WL* sina(l + cosa)coscp 


Following the same procedures as before would yield 


a tt + 


( j2 \ 

J X 
4j 2 

J 


co^ tana = sin a (1 + cos a)] 


(48) 


(49) 


which for small a becomes 


% + 


f J~ x ^ 2 W r . 

4 J 2 J 


* 


= 0 


v 


(50) 


and again, the motion is stable if 


GOjy > 2 


(2 JWL* f 2 

J x 


(51) 


where J is now computed with respect to the support point. 

As it is generally known from the top experiments, a decrease in the angular speed produces an 
increase in the a angle and finally the top falls on the table. It is also known that if the angle a increases 
above a certain limit, the motion destabilizes, an effect that is more important at low speeds. 
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] 

> 

%r 

“3 


L 
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Figure 7. — Bevel gear assembly. 


The same type of analysis may also apply for certain types of bevel gears (fig. 7). When the driver 
gear is lifted, the driven shaft moves around the fixed point (the joint). Here the analysis would somehow 
be different because of the influence of the elastic element. The driven shaft is also subject to the 
influence of a hydrodynamic force due to the fluid film in the joint. 

The study of the motion of bodies with various restraints is an important part of dynamics, and 
different aspects are widely treated in various textbooks on gyroscopes, dynamics of mechanism, flight 
dynamics, rotordynamics, and so on. Closed-form solutions can be found for simplified cases and 
analytical stability considerations can sometimes be deduced. A full review of the analytical solutions is 
way beyond the scope of this paper. However, as an example, we will provide two more results that 
appear in the plain restrained motion of a thin disk. 

Figure 8 illustrates the first example, the rolling of a disk with the disk plane nearly vertical (a « 90°). 
It can be proven (ref. 1) that the disk oscillates around the vertical plane when the spin is larger than 


<4 > ■ 


WJR 
J + mR' 


(52) 


where R is the radius of the disk and m is the mass. As indicated in the appendix, the mass moments of 
inertia are J x ~ mR 2 / 2 and J= mR 2 /4. The stability condition (eq. (52)) is applicable for cylinders where 
the radius is much larger than the height. The classical example of this motion is a coin rolling along a 
table. 

Also, when the disk spins upright (fig. 9), it can be proven (ref. 1) that the spinning motion is stable 
as long as 


i .i 

vf = v > 


WRJ 
J + mR' 


(53) 


The expression of the transverse moment of inertia J has just been indicated above. An example of this 
motion is the classical spinning of a penny or silver dollar on a table top. 
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7 


Video Clip 


The video illustrates the skill of a street merchant walking his empty, open-end, 55 -gal drum about a 
heavily trafficked street with nearly perfect control as one might walk a dog. The drum motion and 
feedback control are maintained by a cadence of regular hand-imparted surface torques at each rotation of 
the barrel. Small angles from vertical require short slaps and large angles require complex upward 
motions of the hand and body. Figure 10 presents some still shots taken in the transient period during the 
lifting of the barrel. A video of this motion can be viewed in the PDF file of this document on the web at 
http://gltrs.grc.nasa.gov/reports/2006/TM-2006-213583/TM-2006-213583.pdf 




(b) 

Figure 10. — (a) Video clip illustrating the dynamics and control of a rotating barrel, 
(b) Lifting the barrel (motion of barrel during transient stage). 


NASA/TM— 2006-213583 


18 




Conclusions 


The dynamic motions of a drum about its center of mass that is restrained by a plane normal to the 
axis of rotation passing through its center of mass at nutation angle a are both entertaining and instructive. 
Yet, once spinning, the motion of the drum can be controlled by regular periodic changes in torque on the 
periphery of the drum, such as once every time the drum rotates. This report presents a preliminary 
analysis of these types of motions. A model is deduced using the theorem of the moment of momentum. 
The governing equations for the nutation angle are deduced. For small values of a, the governing 
equations are simplified and for the illustration, control seems to be an accomplishment of those skilled in 


rhythm and agility. The drum can be stabilized provided that eoj^ > 2 


(. JWL *) 


1/2 


'J 


X 


where eo v is the 


minimum stable speed, J is the transverse mass moment of inertia, W is the weight of the cylinder, L* is 
the distance from the base to the center of mass, and Jx is the mass moment of inertia about the 
longitudinal axis OX of the barrel. 

For the open drum or barrel of the illustration, assuming constant material thicknesses, that rotation 
speed is calculated to be 


co x >16.7 rad/sec with open-drum weight of 143 N and LID = 1.5 (54) 

cox >10.8 rad/sec with open-drum weight of 102 N and LID - 1 (55) 

For the open drum of the illustration, the rotating speed was between 1 and 3 rotations/sec or (6.3 and 
18.8 rad/sec) depending on angle a. 
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Appendix — Moments of Inertia 

Details of Calculations of Mass Moments of Inertia 

The expressions of the inertia moments can be found in various engineering handbooks. Some 
expressions are presented in the body of this work for completeness. The geometric elements are shown 
in figure 1 1 . 

For the plain cylinder or shell, 

Jx s - | |y 2 + z 2 )dm - p| |y 2 + z 2 jdF = pxj |y 2 + z 2 jda ? (56) 

n n s 

where da s = rdBdx, y = r sin 0, and x = r cos 0, so 

2jr L-h 

Jx s = pi | | r 3 d9dx = 2jipxr 3 Z. = InpxrLr 2 - m^r 2 (57) 

0 -h 


Next, 



2nL-h 


+ z‘ 


! )da,=pxri| J 


0 -h 


+ r 2 cos 2 0 ) d 0 d x = 27xpxr 1 


,-h 

( 2 ^\ 

1 

7 7 1 

X + — 
2 

-h 

v 7 


dx 


(58) 


where cos 0 = (cos2 0 +l)/2 and r\ = ( Rod + Rid)/ 2 « R, so 

, n L—h 


Jy s - 27ipxr 1 L 




v 3L y 


+ 


. 2L . i 

k J d-h 


■ mi 


I? -3Lh + 3h 2 
3 


\ 

( ^ 




R~ 



+ 





2 


/ 

v J 

_ 


(59) 



Figure 11. — Geometric elements. 
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For h = L* = qL/2 , where h = L* and t = + m 2 ) are measured from the base, 


Jy s = m x 



( 9 \ 

4-6^ + 3^ 2 

f 

12 

V J A 



( 60 ) 


and for £,= l. 


Jy s = m x 



( 


( 9 "N 



Lr 

+ 

R~ 



12 


2 


_ 

v 2 


v 2 

_ 


(61) 


For end caps or bases, 


R2n 


Jxfr = pi J [y 2 +z 2 j d a/, = pxj J r 3 d r d 0 = 


f 7tpx R 2 R 2 x 


0 0 


( m 2 Rl 1 


(62) 



+ 


jda^, — » pxj 


z 1 da^ 


: px 


R 2n 

II 

0 0 


3 2 

r drcos“ 


OdO = 


f 9 9 \ 

TtpxT? i?“ 


6 ? A 

4 

V 2 


4 

V 2 


(63) 


So, for the plain cylinder, 


Jx s = m | R 


2 


(64) 


and 


Jy s = m x 


L 2 

( 9 \ 

4-6^ + 3 ^ 2 


12 

V 2J 



v j 


and for ^ = 1, the geometric center and mass center coincide: 


Jy s = m x 



A 9 A 


( "> 5 



Zr 

+ 

7?“ 



12 


2 


_ 

V 2 


v 2 

_ 


(65) 


( 66 ) 


For the open drum or cylinder plus the base (or end cap), 


Jx = Jx s + Jxfr = m x R z + 


m 2 R z 


■ m x R z 


1 + m 2 

2m 


1 J 


(67) 


and 
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where, from the moment balance in static equilibrium, = m 2 l{m\ + m 2 ) and \ = m\t(m\ + m 2 ). 
For the drum or barrel or cylinder with the same ends. 



Main Formulations of Mass Moments of Inertia 

For the plain cylinder or for the shell (symmetric), 

{ J Xs)cg =T ^ r2 and ( J ^)cG = ~v~ fir] ^7] (71) 

y z ) v 0 / \ L J 

For the open drum (shell plus end cap, asymmetric), 
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(' J X )cg * “ )cG + (■ Jxb )cG ~ miR ~ 


1 + 


0.57772 
m \ , 


(72) 


and 


(' J )cG * ~( J s)cG +m l 


f Cl?~ 

v 2 , 


( J *)cc +W7 2 


K) 


~l2 


v 4 , 


jr> 2 r - ^ 


+ 2 




+ 


v-w 


777 2 

v m \ j 


R 

\ L J 


+ H-S 




*2 


(73) 


The location of the center of gravity of the open drum ( CG *) from the center of gravity of the shell 
(CG S heii) is defined by 


S =■ 


7?7 2 


( 777J + 777 2 ) 

For the closed drum (shell plus both ends capped, symmetrical), 

{ J x)cg = ( J Xs)cG +2 ( J Xb) C G =m l R 


miR 1 

1 + 

( \ 

m 2 



l m \ )_ 


(74) 


(75) 


and 


(■>,) 


CG 


v 2 , 


(\ 

- 1 + 

v 


r \ 

TTTj 

V m \ J 


fR\ 2 


J 


1 + 


f \ 

TTTj 

v m \ J 


(76) 


The inertial properties of the typical barrel ( L = 0.87 m and D = 0.58 m) are given in table II. 


TABLE II.— INERTIA DATA OF BARREL 
[Length, L, 0,87 m; diameter, D, 0,58 m.] 


Geometry 

Mass moment of inertia about 
longitudinal axis of barrel, J x , 
kg m 2 

Transverse mass moment 
of inertia, J, 
kg m 2 

Weight of cylinder, 
W, 

N 

Drum 

1.225 

2.187 

163.3 

Open drum 

1.138 

1.694 

142.9 

Plain cylinder (shell, no ends) 

1.050 

1.313 

122.5 
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